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ABSTRACT 

We report on the discovery of a new correlation between global parameters of the hot interstellar gas 
in elliptical galaxies. We reanalyze archival Chandra data for 30 normal early-type systems, removing 
the contributions of resolved and unresolved point sources to reveal the X-ray morphology of the hot 
gas. We determine the half-light radius, Rx, and the mean surface brightness, Ix, from the gas 
surface brightness profiles. A spectral analysis determines the temperature, Tx, of the gas within 3 
optical effective radii. We find that the galaxies lie on an X-ray Gas Fundamental Plane (XGFP) 
of the form Tx oc R X 2S I X 22 . This is close to, but distinct from, a simple luminosity-temperature 
relation. The intrinsic width of the XGFP is only 0.07 dex, nearly identical to that of the stellar 
(optical) fundamental plane (SFP). This is surprising since X-ray gas masses are typically ~ 10~ 2 
of the stellar masses. We show that the XGFP is not a simple consequence of the virial theorem or 
hydrostatic equilibrium, and that it is essentially independent of the SFP. The XGFP thus represents 
a genuinely new constraint on the hydrodynamical evolution of elliptical galaxies. 
Subject headings: galaxies: cooling flows — galaxies: elliptical and lenticular, cD — galaxies: ISM — X- 
rays: galaxies — X-rays: ISM 



1. INTRODUCTION 

Elliptical galaxies lie on a two-dimensional locus, 
known as the fundamental plane (FP), in the space 
defined by the optical effective radius R e , optical sur- 
face brightness J R , and stellar velocity disp ersion cjq 
Ipiorgovski fc Davisll987HDressler et al.lll987lh The FP 
is understood to be a consequence of the virial theo- 
rem, modified by stellar population variations and struc- 
tural nonhomology for systems of different luminosity 
ijTruiillo et all I2004L and references therein). FP-like 
scaling relations have been found for oth er types of sys- 
tems as well, including galaxy clusters (|Schaeffer et alJ 
1993). Several authors have pointed out "cluster funda- 
mental pla ne" relations involv ing either optical parame- 
ters alone llAdami e t al.||19 98[) or mixtures of X-ray and 
optical parameters (Fritsch & Buchert 1999; Miller et al. 
1999). Because the baryonic mass in clusters is dom- 
inated by intracluster gas at temperatures > 1 keV, X- 
rays are a better tracer of o rdinary matter than starlight. 
iFuiita fc Takaharal l)1999ajl show the existe nce of a clus - 
ter FP purely in X-ray parameters (see also lAnnislll994f) . 
which th ey show can be consistent with simple spherical 
collapse l|Fuiita fc Takaharalll999b]l . 

An "X-ray fundamental plane" for elliptical galaxies — 
actually a relation between stellar dispersion, X-ray lu- 
minosity, and half-l i ght ra dius — has been suggested by 
iFukueita fc Peeblesl (^)99). But in normal ellipticals, 1 



the mass of X-ra y gas is typically only a few percent of 
the stellar mass ijBreeman et al.lll992|) . Much of the X- 
ray flux comes from low- mass X-ray binaries (LMXBs); 
if not accounted for, this guarantees a correlation be- 
tween X-ray and optical luminosity. Normal ellipticals 
are tho ught to lose most o f their gas in supernova-driven 
winds ijCiotti et al.lll99l|) . and high-resolution observa- 
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1 For our purposes, normal ellipticals are those that are not at 
the centers of cluster potential wells. 



tions using Chandra show that many appear disturbed, 
suggesting redistribution of gas through shocks, nuclear 
activity, merger s, or interaction with the intergalac- 
tic medium (e.g iFinoguenov fc Jonesl 120011: IJones et alJ 
l2CMlMachacek et alJl2005HStatler fc McNamarall2002|) . 
Thus, the visible hot gas is only a tenuous leftover of a 
complex hydrodynamical history, and may be far from 
equilibrium. 

Nonetheless, in this Letter we report the discovery of a 
fundamental plane relation for the X-ray gas alone. We 
show that the X-ray Gas Fundamental Plane (XGFP) is 
distinct from — but as tight as — the stellar fundamental 
plane (SFP). Unlike the SFP, the XGFP is not a simple 
consequence of the virial theorem. In fact, the XGFP 
is almost completely decoupled from the SFP, and thus 
constitutes a new constraint on the evolution of normal 
elliptical galaxies. 

2. DATA REDUCTION AND ANALYSIS 

2.1. Chandra Archive Sample and Pipeline Reduction 

We analyze a sample consisting of 56 E and E/SO 
galaxies having non-grating ACIS-S exposures longer 
than 10 ks in the Chandra public archive. Brightest clus- 
ter galaxies and objects with AGN-dominated emission 
are excluded. All observations are uniformly reprocessed 
using version 3.1 of the CIAO software and version 2.28 
of the calibration database. Flares are removed by itera- 
tively applying a 2.5cr threshold. For the quiescent back- 
ground, intervals more than 20% above the mean count 
rate are excised, to match the blank sky background 
fields. We restrict photon energies to the range 0.3- 
5keV, further divided into soft (0.3-1.2 keV) and hard 
(1.2-5.0 keV) bands. Monoenergetic exposure maps are 
created in steps of 7 in PI (~ 100 eV). An image is ex- 
tracted for each 14.6eV-wide PI channel, and divided 
by the energetically closest exposure map to create a 
photon-flux-calibrated "slice." The slices are summed to 
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produce the calibrated photon flux images. Full details 
of the analysis and results, alo ng with data products, wi ll 
be presented in future papers ijDiehl fc Statlerll2006alblh 

2.2. Isolating the Gas Emission 

Point sources are identified in each band by CI AO's 
wavdetect tool. Regions enclosing 95% of the source 
flux are removed and the holes filled with counts drawn 
from a Poisson distribution, whose pixel-by-pixel expec- 
tation value is determined by adaptive interpolation us- 
ing the asmooth tool in the XMMSAS package. A uni- 
form background is determined from fits to the radial 
surface brightness profile, and subtracted. 

To remove the contribution of unresolved point 
sources, we use the fact that the hot gas and LMXBs 
contribute differently to the soft and the hard bands. 
Let S and H represent the background-subtracted soft 
and hard images. We can express both in terms of the 
unresolved point source emission P, the gas emission G, 
and their respective softness ratios 7 and 5: 



S = jP + SG 
H=(l-j)P + (l-S)G. 

The uncontaminated gas image is then given by 
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Assuming that resolved and unresolved LMXBs share 
spatially independent spectral properties, we can use the 
resolved sources to determine the constant 7. We take 
sources between 5" and 5R e from the center, excluding 
high luminosity sources (> 200 counts). For systems with 
> 10 sources meeting these criteria, we fit an absorbed 
power law model to the combined point source spectrum, 
with hydrogen column density fixed at the Galactic value 
for the line of sight. Integrating the model over the soft 
and hard bands yields 7. For other galaxies , we exploit 
the u niversal nature of the LMXB spectrum ijlrwin et alJ 
2003) . A simultaneous power-law fit to all low-luminosity 
(Lx < 5 x 10 37 erg s" 1 ) point sources in our sample gives 
a photon index of 1.603. This model is used to derive 7 
for the source-poor galaxies. 

The coefficient S is determined si milarly, from the fit 
of an APEC thermal plasma model l|Smith et alJl20Q]Tl 
to the hot gas emission (see § I2.3[) . This approach as- 
sumes isothermal gas throughout the galaxy. In case of 
a temperature gradient, one would have to account for 
the spatial dependence of 5, but this approach is beyond 
the scope of this Letter. 

2.3. Physical Parameters for the Gas 

To produce a radial profile, we adaptively bin the gas 
image G into circular annuli. In 12 cases there is insuffi- 
cient signal to fit the spatial profile. We fit the remain- 
ing profiles with Sersic models to derive X-ray half-light 
radii Rx and mean enclosed surface brightnesses Ix 2 
We discard 14 objects with Rx larger than the size of 
the observed field. The final sample of 30 objects has B 
absolute magnitudes in the range —22.5 < Mb < —19. 
Twenty-one are group members, and 14 are brightest 



group galaxies l)Garcialli993|) . Eight of the 14, plus 3 
additional objects, are centr al members of X-ray- bright 
groups in the GEMS survey ijOsmond fc Ponmanl l2004) . 
We find that up to 55% (typically 10-30%) of the dif- 
fuse emission in the final sample comes from unresolved 
LMXBs. 

The spectrum of the diffuse emission between 0.3 and 
5 keV is extracted from a circular region 3i? e in radius, 
excluding resolved point source regions. We fit the spec- 
trum using the SHERPA package, adopting a single tem- 
perature APEC thermal plasma model for the hot gas 
and a power law for the unresolved LMXBs. The nor- 
malizations of both components, the gas temperature 
Tx, and (in most cases) gas metallicity are allowed to 
vary. The redshift is set to the value give n in the Lyon- 
Meud on Extragalactic Database (LEDA; iPaturel et alJ 
Il997|) . For low signal-to-noise spectra, the abundances 
are held fixed at the solar value. The photon index of 
the power law component is determined by a simultane- 
ous fit to the spectrum of the lowest luminosity resolved 
point sources (§ I2.2fl . Single-temperature APEC models 
are poor fits (reduced y 2 > 2) to 12 objects. In these 
cases the temperature should be interpreted as an emis- 
sion weighted average. Excluding these galaxies from the 
sample does not change the results. 

All errors are assumed to be Gaussian, described by 
a covariance matrix. Statistical errors on Rx, Ix, and 
Tx are obtained from the fitted models. We adopt dis- 
tances and errors obtained from surface brightness fluctu- 
ations l)Tonrv et al]l200ll) where available; otherwise we 
use LEDA values 3 and assume errors of 15%. We adopt 
a 10% uncertainty in 7 from the galaxy-to-galaxy scatter 
in the photon indices fitted to the composite point-source 
spectra, and a 5% uncertainty in S from a comparison of 
values obtained from the spectral fits and from direct 
integration of gas-dominated spectra over the hard and 
soft bands. We measure the effect on Rx and Ix by 
repeating the spatial fits with altered values of 7 and S. 

Tabled lists the base-10 logarithms of Rx (kpc), Ix 
( erg s" 1 cm~ 2 arcsec -2 ) and Tx ( keV) for the final sam- 
ple of 30 galaxies, and the corresponding non-zero ele- 
ments of the covariance matrix. 

3. THE X-RAY GAS FUNDAMENTAL PLANE 

The distribution of galaxies in the parameter space 
(\ogRx,\ogIx,\ogTx) is nearly planar. This is seen 
clearly in Figure ^ which shows face-on and edge-on 
views of the XGFP. We determine the intrinsic distribu- 
tion in this space by fitting a probability density in the 
form of a tilted slab with finite Gaussian width, taking 
correlated errors into account. 

We can express the XGFP in the form 



T x cx R x I b x , 



(4) 



where a and b determine the orientation of the plane. 
The best-fit values are a = 0.28 and b = 0.22 (Figure 
|2| • The formal 1-dimensional errors in a and b are 0.045 
and 0.037, respectively; but as the figure shows, the er- 
rors are correlated. The intrinsic width of the XGFP 
is very small, with a value of .068 ± 0.012 dex, wh ich 
is identical to that of the SFP l|Bernardi et alJl2003l) to 



2 Fits using /3 and double-/? models are unphysical in more than 
half the cases, implying divergent fluxes at large radii. 



3 LEDA distances are obtained from a B-band Faber-Jackson 
relation of the form Mg T = — 6.21og<r — 5.9. 
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TABLE 1 

Physical parameters for the X-ray gas 



— 

ID 




log 




log 




log 


— 

Crr 




Cu 


— 

Ctt 


— — 

Cm 




Rx 


Ix 


T x 


xlO 4 


xlO 4 


xlO 6 


x 10 4 


11459 


0.93 


— 16.73 


—0.318 


276.3 


48.4 


259.5 


—29.2 


14296 


0.35 


— 15.18 


—0.055 


59.3 


32.8 


56.2 


— 12.1 


in lyo 


1 ^ZL 
1 .o*± 


1 fi 7Q 


—0 114 


44 9 


zu.u 


'iO.O 


— u.u 


N315 


0.62 


— 15.77 


—0.196 


60.5 


58.1 


75.4 


— VIA 


N533 


0.73 


— 14.99 


—0.009 


53.1 


90.0 


11.4 


— 13.6 


N720 


0.96 


— 16.46 


—0.247 


48.0 


41.3 


36.9 


—20.0 


N741 


0.64 


— 15.27 


—0.016 


84.0 


91.5 


72.8 


—30.5 


N1 404 


0.52 


— 15.31 


—0.234 


17.5 


40.7 


7.6 


—1.6 


N1407 


0.56 


— 15.61 


—0.061 


29.8 


32.5 


6 4 


— 1.4 


N1553 


1.02 


—17.10 


—0.392 


88.2 


1236.5 


190.2 


2.1 


N2434 


1.16 


— 17.42 


—0.273 


362.9 


29.7 


621.3 


— 79.1 


N3923 


0.23 


— 15.32 


—0.322 


47.0 


24.1 


217.4 


—7.2 


N4125 


1.09 


— 16.82 


—0.356 


61.2 


38.0 


57.6 


—6.2 


N4261 


0.22 


— 15.27 


—0.110 


21.0 


46.2 


56.4 


—5.4 


N4374 


0.28 


-15.30 


-0.151 


6.1 


65.0 


15.4 


-0.9 


N4526 


0.58 


-16.71 


-0.450 


405.6 


95.5 


1328.8 


-89.2 


N4552 


0.01 


-15.20 


-0.245 


9.0 


545.0 


20.1 


-2.2 


N4621 


1.29 


-18.05 


-0.629 


456.5 


74.8 


2359.7 


-133.9 


N4636 


0.76 


-15.60 


-0.160 


73.1 


1092.3 


0.9 


-98.1 


N4649 


0.35 


-15.08 


-0.096 


15.8 


39.4 


2.0 


-7.3 


N4783 


1.14 


-16.57 


0.053 


479.8 


114.8 


878.3 


-105.1 


N5044 


1.07 


-15.24 


-0.041 


34.1 


55.8 


1.5 


-3.1 


N5846 


0.91 


-15.74 


-0.152 


17.3 


40.0 


4.3 


-0.5 


N6482 


1.54 


-16.34 


-0.131 


286.8 


118.1 


21.3 


-67.3 


N7052 


1.08 


-16.13 


-0.278 


1323.9 


507.3 


626.0 


-209.1 


N7618 


1.59 


-16.34 


-0.095 


73.8 


27.2 


35.4 


-6.8 




- 1 7 -IB -15 

log( R x - a:is h "■ 95 T x »•« ) 

Fig. 1. — Face-on (top) and edge-on (bottom) views of the X-ray 
Gas Fundamental Plane. Axis labels indicate eigenvector compo- 
nents. Symbol sizes roughly indicate the relative positions into and 
out of the page. Error bars indicate lcr-projections of the covari- 
ance matrices. Arrows illustrate the sense of view relative to the 
fundamental measured parameters. 



within the errors. The fit is robust and not sensitive to 
the choice of model for the surface brightness profile. Ex- 
cluding the brightest group galaxies or galaxies with bad 
single-temperature fits results in a and b values within 
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Fig. 2. — Exponents a and b, describing the orientation of 
the XGFP according to equation I4"l . Cross and ellipses indicate 
best-fit values and confidence regions, respectively. Diagonal line 
marks combinations of a and b corresponding to pure luminosity- 
temperature correlations of the form L oc T" . 



the 68.3% confidence ellipse in Figure |21 

4. DISCUSSION 
4.1. Relation to Known Scaling Laws 

If a and b obeyed the relation a = 2b, the XGFP would 
be equivalent to a simple luminosity-temperature rela- 
tion of the form Ljf.gas oc T x - , indicated by the solid line 
in Figure The model fit rules out a pure Lx, gas -Tx 
relation at > 99.7% confidence. The Lx, gas -Tx relation 
represents a nearly edge-on vi ew of the XGFP, analog ous 
to the Faber-Jackson relation ijFaber fc ,Iacksonll97 6') for 
the SFP. As in the optical case, the XGFP accounts for 
much of the intrinsic scatter (0.091 dex) in the Lx,g& s -Tx 
relation. 4 For a given Tx, galaxies with more extended 
gas emission are more luminous than compact objects. 

The Lx~Tx relation most nearly consistent with the 
XGFP is given by £x, g as oc T 8 5 (Figure^. A naive prin- 
cipal component analysis in t he L x . ? as -Tx plane yields 
a shallower exponent of 5.9. Q : Sulliv an et alJ (2003) 
find Lx.tot oc T x ' s±0 ' 7 for their complete sample, and a 
steeper exponent of 5.9 ± 1.3 when they exclude galaxies 
with prominent temperature gradients. However, their 
X-ray luminosities include unresolved LMXBs, the re- 
moval of which would steepen the Lx~Tx relation. 

A compariso n can also be made with th e Lj-a relation 
for ellipticals. iMahdavi & Gellerl (|2001f) find L x ,tot oc 



<7 - 1 - 6 , and predict this relation to steepen to Lx. ; 



oc 



r 12±5 



if Lx is restricted solely to the hot gas. Us- 
ing the temperature-d ispersion correlation a oc T x 56 
ijO'Sullivan et al.ll2003D . we can approximate our clos- 
est Lx~Tx relation as £x,gas oc a 15 , which is consistent 

w ith the earlier resul t. 

iFuiita fc Takaharal l)1999a|) obtain a result vaguely sim- 
ilar to ours for clusters of galaxies. They find an X-ray 
cluster FP connecting core radius, central density, and 
mean cluster temperature. Assuming a constant value 
of f3 — 2/3 for their surface brightness profile fits, we 
can translate their cluster FP to our parameters, finding 
Their relation is significantly inclined 



T; 



A' 



oc R X 57 I° X 32 



4 Gas luminosities are obtained by summing the observed flux 
in the field of view and using the fitted Sersic law to account for 
the flux at larger radii. 



4 



Diehl & Statler 



to our XGFP, and is close to the relation L a T 3 . How- 
ever, the cluster FP deviates from a pure L-T relation 
in a manner similar to that of the XGFP. 

iFukugita fc Peebles! l)1999l) suggest a galaxy fundamen- 
tal plane with mixed X-ray and optical parameters, Lx, 
Rx and a 2 . Their sample consists of 11 galaxies with 
ROSAT, Einstein or A SCA observations. We do not 
reproduce their result with our larger Chandra sample. 
The reason for this discrepancy may be that their X-ray 
luminosities are corrected neither for light outside the 
field of view nor for the contribution of point sources. 

4.2. Independence of the XGFP and the SFP 

The virial theorem connects a system's total mass M 
with its characteristic radius Rm and dispersion o~m- 
This relation produces an observable SFP because mass 
maps to luminosity by way of the mass-to-light ratio, 
and the stellar R e and a are surrogates for Rm and ctm- 
For the XGFP to be another manifestation of the virial 
theorem, one would require similar mappings from mass 
parameters into X-ray observables. One example is hy- 
drostatic equilibrium, which links Rx and Tx for a given 
potential. Others might connect gas mass to total mass 
or to a measure of gas retainability such as a 2 /Tx- How- 
ever, these relations would be identifiable in correlations 
between optical an d X-ray parameters. O nly the known 
Tx~<y relation (e.g. lO'Sullivan et al.ll2003ft and an addi- 
tional, very weak Tx~R e relation are supported by the 
data. X-ray gas masses are uncorrelated with Mb and 
other optical properties. 

If the SFP and XGFP were linked, then each plane 
would represent a projection of a higher dimensional, 
more fundamental hyperplane into the corresponding 3- 
parameter subspace. We test this hypothesis by ana- 
lyzing the (logRx, log/x, logTx, logi? e , log/ e , logo- 2 ) 
space with principal component analysis (PC A). In this 
space our sample is reduced to 25 objects with reliable 
X-ray and optical parameters. We define the 3-vectors 
nx and no to be the normals to the fundamental planes 
in X-ray and optical parameters, respect ively. From the 
SFP we have n = [0.69,0.51,-0.51] (|Bernardi et all 
1200^ and from the XGFP, n x = [0.26,0.21,-0.94]. 
If the planes are completely independent, the eigenvec- 
tors corresponding to the two smallest eigenvalues will 
have the form Ni — [/i nx, y/l — M 2 no] and N 2 = 



[a/1 — fi 2 nx, — M no]- Here, N± and N 2 are 6- vectors, 
and /j, can be any number between —1 and 1, depend- 
ing on the relative scatter in the two planes. We obtain 
Nx = [-0.21,-0.16,0.84,0.20,0.15,-0.39] and N 2 = 
[0.01, 0.03, -0.38, 0.79, 0.38, -0.28] . This result is not far 
from the above prediction with /i = —0.40. Furthermore, 
the scatter about both planes is not reduced by going to 
higher dimensions. Both results point to, at most, weak 
coupling of the subspaces. 

Optical and X-ray parameters are known to be coupled 
through the observed Tx~cr 2 relation. We test whether 
this coupling could have a measurable effect on the eigen- 
vectors, using Monte Carlo simulations of 25-object sam- 
ples following the fitted SFP and XGFP relations. We 
take the SFP and XGFP to be independent, except for 
an intrinsic linear correlation between log Tx and log a 2 
with a variable amount of scatter. We consider three 
coupling strengths: none, weak, and strong. Weak cou- 
pling reproduces the observed width (0.10 dex) of the 
Tx-cr 2 relation and, by way of the SFP, the somewhat 
weaker Tx~R e relation; strong coupling reduces the Tx- 
cr 2 width by a factor of 3. We perform PCA on the sim- 
ulated 6-d data and measure the alignment of the best 2 
eigenvectors with those obtained from the real data. All 
3 cases reproduce the eigenvectors to within the Poisson 
noise, showing that the observed Tx-cr 2 relation does not 
change the relation between a decoupled XGFP and SFP 
at a level that can be resolved with a 25-object sample. 

The data are thus consistent with the XGFP and SFP 
being almost completely independent. The XGFP can- 
not be understood as a simple consequence of the virial 
theorem or hydrostatic equilibrium. Instead, the XGFP 
represents a new constraint on the hydrodynamic evolu- 
tion of elliptical galaxies. 
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